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Abstract 

In this paper we investigate for further symmetry properties of the nonlinear fin equations of the general form ut = (E(u)u x )x + 
h(x)u rather than recent works on these equations. At first, we study the projective (fiber-preserving) symmetry to show that 
equations of the above class can not be reduced to linear equations. Then we determine an equivalence classification which 
admits an extension by one dimension of the principal Lie algebra of the equation. The invariant solutions of equivalence 
^"*^ , transformations and classification of nonlinear fin equations among with additional operators are also given. 
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■ 1 Introduction 

in 

t***'- ■ Investigations for symmetry properties of mathematical models of heat conductivity and diffusion processes [9] 
C*~) ' are traditionally formulated in terms of nonlinear differential equations which often envisage us with difficulties in 
qq , studying. To solve this problem, symmetry methods play a key role for finding their exact solutions, similar solutions 
f*^; ■ [2,3,15,17,18] and invariants. 
G\ . 

In this study we generalize the study of a class of nonlinear fin equations which has been recently studied in some 
• • references and specially in [17,18]. So, we are dealing with the class of nonlinear fin equations of the general form 
^ ■ 

u t = (E(u)u x ) x + h(x)u, (1.1) 

0$ ' in which we assumed that E u ^ 0, u is treated as the dimensionless temperature, t and x the dimensionless time and 
space variables, E the thermal conductivity, h = —N 2 f(x), N the fin parameter and / the heat transfer coefficient 
[2]- 

The Lie point symmetry in linear and nonlinear case, the condition E u — corresponds to the linear case, 
the class of nonlinear one-dimensional diffusion equations when h = 0, the class of diffusion-reaction equations 
when h = cons., the case which the thermal conductivity is a power function of the temperature and additional 
equivalence transformations, conditional equivalence groups and nonclassical symmetries have all investigated and 
listed in Table 1 of [18]. The point symmetry group of nonlinear fin equations of class (1.1) were considered in a 
number of papers (e.g. see [2] for the physical meaning and applications of the equation). The Lie algebra of the 
point symmetry group of Eq. (1.1) is 

* -(!), d.2) 
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In the next section, we concern with the problem of finding projective symmetry group of Eq. (1.1) as a special case 
of the point symmetry group; since it may have important information about the equation. 

The equivalence classification of Eq. (1.1) in the special case 
u t = (E(u) u x ) x , (1.3) 
has performed by L.V. Ovsiannikov [14]. In [17] authors carried out the another special case of equivalence group 

t = 5 1 t + S 2 , x = S 3 x + S 4 , u = S 5 u, E = 5^ 1 5lE, h = 5^ 1 h, (1.4) 

of Eq. (1.1) when Si, i = 1, • • • , 5, are arbitrary constants and Si S3 S5 ^ 0. The more general class of nonlinear fin 
equations is the nonlinear heat conductivity equations of the form 

u t = F(t, x, u, u x )u xx + G(t, x, u, u x ), (1.5) 

which admits non-trivial symmetry group. The group classification of (1.5) is presented in some references [1,11] . 

However, since the equivalence group of (1.5) is essentially wider than those for particular cases, the results of [1,11] 
cannot be directly used for symmetry classification of particular ones. Nevertheless, these results are useful for finding 
additional equivalence transformations in the class of our problem. Therefore in contrast to the above works, in the 
last two sections of this paper, we study group classification of Eq. (1.1) under equivalence transformations in the 
general case. Furthermore, a number of nonlinear invariant models which have nontrivial invariance algebras are 
obtain. 

From [19] we know that if the partial differential equation possesses non-trivial symmetry, then it is invariant 
under some finite-dimensional Lie algebra of differential operators which is completely determined by its structural 
constants. In the event that the maximal algebra of invariance is infinite-dimensional, then it contains, as a rule, 
some finite-dimensional Lie algebra. Also, if there are local non-singular changes of variables which transform a 
given differential equation into another, then the finite-dimensional Lie algebra of invariance of these equations 
are isomorphic, and in the group-theoretic analysis of differential equations such equations arc considered to be 
equivalent. To realize the group classification, we use of the proposed approach consists in the implementation of an 
algorithm explained and performed in references [1,10,14,16]. For this goal, our method is completely similar to the 
way of [6] for the nonlinear wave equation u u — f(x, u)u xx + g(x, u). 



2 Projective symmetries of nonlinear fin equations 

In this section, we are concerning with group classification of nonlinear fin equations by projective transformations 
group as a special case of the point symmetry group. Our study is based on the method of [13] for Lie infinitesimal 
method. 

The equation is a relation among with the variables of 2-jet space J 2 (R 2 ,R) with (local) coordinate 

(t,x,u,u t ,u x ,u tt ,u tx ,u xx ), (2.6) 

where this coordinate involving independent variables t, x and dependent variable u and derivatives of u in respect to 
t and x up to order 2 (each index will indicate the derivation with respect to it, unless we specially state otherwise). 
Let M be the total space of independent and dependent variables resp. t, x and u. The solution space of Eq. (1.1), (if it 
exists) is a subvariety 5a C J 2 (R 2 ,R) of the second order jet bundle of 2-dimensional sub-manifolds of M. If we wish 
to preserve the bundle structure of the space A4, we must restrict to the class of fiber-preserving transformations in 
which the changes in the independent variable are unaffected by the dependent variable. Projective or fiber-preserving 
symmetry group on M is introduced by transformations in the form of 

x = (j){t,x), t = x(t,x), u = ij)(t, x,u), (2.7) 

for arbitrary smooth functions 4>, %, ip. Also assume that 

v := e (t, x)^ + fit, x)^+ V(t, x > u )^ ( 2 - 8 ) 
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with coefficients as arbitrary smooth functions be the general form of infinitesimal generators which signify the Lie 
algebra q of the projective symmetry group G of Eq. (1.1). The second order prolongation of v [13,14] as a vector 
field on J 2 (R 2 ,R) is as follows 

v {2) ■■=v + r 1 t ^-+v x ^-+V tt ^-+V tx 7 ^-+V xx ^-, (2-9) 
ou t ou x ou tt ou tx ou xx 

where rf ,rf and rf*, n tx , rf x are arbitrary smooth functions depend to variables t,x,u,p,q and (2.6) resp. These 
coefficients are computed as following 

V J =Vj(Q) + Z 1 u Jit + Z 2 uj, x (2.10) 

where T> is total derivative, J is a multi-index with length 1 < | J| < 2 of variables t, x and Q = u — £} u t — £ 2 u x 
is characteristic of v [13]. According to [13], v is a projective infinitesimal generator of Eq. (1.1) if and only if 
«( 2 )[Eq. (1.1)] = 0. By applying on the equation we have the following equation 

£ 2 h x u + r](D uu q 2 + D u q x ) - rf + 2rf D u q + r\ xx D = 0, whenever Eq. (1.1) is satisfied. (2.11) 

In the extended form of the latter equation, functions £ , £ 2 and r\ only depend to t, x, u rather than other variables, 
i.e. u t ,u x , u t t,u tx , u xx , hence the equation will be satisfied if and only if the individual coefficients of the powers of 
Ut, u x and their multiplications vanish. This tends to the following over-determined system of determining equations 

e u £ = o, e£ = o, r) u - e t +e£ x = o, 

E u r 1 + Er ]u -2Eg=0, £ 2 h x u - r) t + Erj xx + hr) = 0, (2.12) 

E uuV + 2 E u ( Vu -&) + E ^ uu = 0, i 2 t +2E uVx + E(2 Vxu - £ x ) = 0. 

By solving Eq. (2.12), the general solution to these differential equations for ^,^ 2 and 77 will be found: 

£{t,x) = c, e{t,x)=0, V (t,x,u)=0, (2.13) 



with arbitrary constant c. Therefore the Lie algebra g spanned by projective infinitesimal generators of Eq. (1.1) is 
id_\ 



Q = {■§;} and the projective symmetry group is nothing but the time translation group. 



Theorem 1 A complete set of all infinitesimal generators of nonlinear fin equation up to projective transformations 
admits the structure of one-dimensional Lie algebra and so is isomorphic to R. The point symmetry group and 
projective symmetry of nonlinear fin equation are equal. 

It is well-known that the existence of a non-fiber-preserving symmetry usually indicates that one can significantly 
simplify the equation by some kind of hodograph-like transformation interchanging the independent and dependent 
variables. But in the case of our problem, we can not use this advantage for simplifying nonlinear fin equations. 

According to [13], any system of partial differential equations which has only a finite-dimensional symmetry 
group is certainly not linearizable, that is, for every change of variables, it can not be mapped to an inhomogeneous 
form of the linear system T>[u] = f, where V is & second order linear differential operator, u indicates dependent 
variables and / denotes smooth functions of independent variables. 

Remark 2 A system of nonlinear fin equations in the form of Eq. (1.1) can not be reduced into an inhomogeneous 
form of a linear system. 



3 Equivalence transformations 

At the present section we follow the method of L.V. Ovsiannikov [14] for partial differential equations. His approach 
is based on the concept of an equivalence group, which is a Lie transformation group acting in the extended space 
of independent variables, functions and their derivatives, and preserving the class of partial differential equations 
under study. It is possible to modify Lie's algorithm in order to make it applicable for the computation of this group 
[1,10,14,16]. Next we construct the optimal system of subgroups of the equivalence group. 
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An equivalence transformation is a non-degenerate change of the variables t, x, u taking any equation of the form 
(1.1) into an equation of the same form, generally speaking, with different E(u) and h(x). The set of all equivalence 
transformations forms an equivalence group G. We shall find a continuous subgroup Gc of it making use of the 
infinitesimal method. 

We investigate for an operator of the group Gc in the general form 

Y := ^(t, x) -| + £ 2 (t, x) A + r)( t , x,u)^- + V (t, x, u, £?, fc) A + X (t, x, u, E,h)^-. (3.14) 

from the invariance conditions of Eq. (1.1) written as the system 

u t = (E(u) u x ) x + h(x) u, E t = E x = 0, h t = h u = 0, (3.15) 

where we assumed that u, E, h arc differential variables: u on the base space (t, x) and E, h on the total space (t, x, u). 
Also in Eq. (3.14) the coefficients depend to variables t, x, u and the two last ones, in addition, depend to E, h. The 
invariance conditions of the system (3.15) are 

Y [«* - (E(u) u x ) x - h{x) u] = 0, Y [E t ] = Y [E x ] = 0, Y [h t ] = Y [h u ] = 0, (3.16) 
where 

Y := Y + 71* — + rf — + 7i u — + rf x — + „ xx -A_ + ^ — + <p x — + v* — + v u — . (3.17) 

is the prolongation of the operator (3.14). Coefficients rj J for multi-index J (with length 1 < |J| < 2) have given 
in section 2 and by applying the prolongation procedure to differential variables E, h with independent variables 
(t, x,u) we have 

<p* = Vtitp) - Et Vt(e) - E x V t {e) - E u Vtin) = V t (<p) - E u V t (j,), 

v x = V x {tp) - E t V x {^) - E x y x (?) - E^V X ( V ) = V x (ip) -E U V X { V ), (3 lg) 

x* = Vt(x) - h t Vt{£) - h x v t (e) - h u v t {r,) = v t {x) - h x v t (f), 
x u = v u {x) -h t v u {e)-h x v u (e) - K v u {n) = v u ( x ) - h x vjf), 

where in view of Eq. (3.15) we have 

~ Odd d 
Vt : =di +Et dE +ht dh = di> 

V x := — +E x — +h x ^ = iz + h xiz' ( 3 - 19 ) 
~ d d 
Vu '-^ + Eu dE^ 

Substituting (3.17) in (3.16) we tend to the following system 

ipu^ + hrj + xu-rf + 2E u u x r] x + tp u u 2 x + Er 1 xx = 0, (3.20) 
Lp l = 0, y x = 0, (3.21) 

X* = 0, X u = 0. (3.22) 



d 


d 


. d 


dh ~ 


dx 




d 


d 










dh ~ 


du 


^ Eu dE 



Replacing relations r\ J (for multi index J with length 1 < \J\ < 2) and x l iX x m Eqs. (3.20)-(3.22) and then 
introducing the relation u t — {E{u) u x ) x + hu to eliminate u t , we have five relations which are called determining 
equations. The four last ones are the determining equations associated with Eqs. (3.21), (3.22), i.e., 

<Pt - E u T) t = 0, ifx + h x tp h - E u r] x = 0, xt ~ h x £ =0, Xu + E u xe ~h x (,l = 0. (3.23) 
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But these relations must hold for every E and h and this fact results in the following conditions 

& 2 = £ = 0, Vt = Vx = 0, (pt = tp x = ip h = 0, xt = Xu = Xe = 0, (3.24) 

so, we find that 

e = S, 2 (x), V = v(u), <p = <p(u,E), X = x(x,h). (3.25) 

Adding these conditions to the remained determining equation, since u t ,u x ,utt,u tx ,u xx are considered to be inde- 
pendent variables, we lead to the following system of equations 

hrj + X u + hu - r] u ) = 0, E u { Vu + & - 2^) +E Vuu + E uu r? = 0, 

+ + -2g) = 0, e£ = o, e£ = o, e£ x = o. (3 - 26) 

This system follows 

£ 1 = 2cit + c 2 , Z 2 = Cl x + c 3 , n = 0, <p = e- u F(E), X = -2c 1 h, (3.27) 

with arbitrary function F = F(E) and constants Ci, 02,03. Therefore the class of Eqs. (1.1) has an infinite continuous 
group of equivalence transformations generated by infinitesimal operators 

* = w Y * = 2 4 +x i- 2h l> y ^^ f ^je- ^ 

Moreover, in the group of equivalence transformations are included also discrete transformations, i.e., reflections 
ti — > -t, xi — ► -x, iti — y -u, Ei — y-E, hi — ► -h. (3.29) 

The communication relations between these vector fields is given in Table 1. The Lie algebra g := {Yp,Yi : i = 
1, • • • , 3 ) is solvable since the descending sequence of derived subalgebras of q: q D = { 2 Yi, Y<i ) D g' 2 -* = {0}, 
terminates with a null ideal. But its Killing form: K(v, w) — tr(ad(i>) o ad(u>)) = 5 03 63 for each v — Vifi + v^Yp 
and w = J2j w jYj + w 4 Yp in g is degenerate and hence g is neither semisimple nor simple. 

Theorem 3 Let Gi be the one-parameter groups generated by the Yi, then we have 

G\ : (t, x, u, E, h) 1 — y (t + s,x,u,E,h), Gi : (t, x, u, E, h) 1 — ► (t,x + s,u, E,h), 

( 0.0U J 

G 3 : (t,x,u,E,h) 1 — ► (te 2s ,xe s ,u,E,he 2s ), G4 : (t,x,u,E,h) 1 — ► (t,x,u,E,h), 

when E is the solution of equation da/F(a) — e'"s for a constant c. Furthermore, if u — f(t,x) for functions 
E and h be a solution of nonlinear fin equation, so are the following functions 

ui = f(t + s,x,u), u 2 — f(t,x + s,u), (3.31) 

for the same functions E and h, u 3 = f(t e 2s , xe s ,u) for the same functions E and h = he 2s and also U4 = f(t, x, u) 
for E and the same h. 

4 Preliminary group classification 

In many applications of group analysis, most of extensions of the principal Lie algebra admitted by the equation 
under consideration are taken from the equivalence algebra Q£. These extensions are called £ -extensions of the 
principal Lie algebra. The classification of all nonequivalent equations (with respect to a given equivalence group 
Gs) admitting Z -extensions of the principal Lie algebra is called a preliminary group classification [10]. We can take 
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Table 1 

Commutator table 



M Yi Y2 n n 

Yi 2Fi 

Y 2 Y 2 

y 3 -2Y1 -y 2 

y 4 



Table 2 
Adjoint table 



Ad 


Fi 


Y 2 


Y 3 


Yi 


Fi 


Y 1 


Y 2 


Y 3 -2sY 1 


Y 4 


Y 2 


Y 1 


Y 2 


Y 3 -sY 2 


Y 4 


Y 3 


e 2s n 


e s Y 2 


Y 3 


Y 4 


Y 4 


Y l 


Y 2 


Y 3 


Y 4 



any finite dimensional subalgcbra (desirable as large as possible) of an infinite-dimensional algebra with basis (3.28) 
and use it for a preliminary group classification. We select the subalgebra jj 4 spanned on the following operators: 

d d d d d d 



It is well-known that the problem of classifying invariant solutions is equivalent to the problem of classifying sub- 
groups of the full symmetry group under conjugation in which itself is equivalent to determining all conjugate 
subalgebras [13,14]. The latter problem, tends to determine a list (that is called an optimal system) of conjugacy 
inequivalent subalgebras with the property that any other subalgebra is equivalent to a unique member of the list 
under some element of the adjoint representation i.e. f) Ad(<?) f) for some g of a considered Lie group. Thus we will 
deal with the construction of the optimal system of subalgebras of q 4 . 



The adjoint action is given by the Lie series 



Ad(cxp( s Y t ))Y 3 =Yj-a [Yi, Y 3 ] + - [Y u [Y h Yj] 



(4.33) 



where s is a parameter and i,j — 1, • • • ,4. The adjoint representations of 04 is listed in Tables 2; it consists the 
separate adjoint actions of each element of g 4 on all other elements. 

Theorem 4 An optimal system of one- dimensional Lie subalgebras of nonlinear fin equation (1.1) is provided by 
those generated by 



1) A 1 =Yi= d t , 

2) A 2 = Y 2 = d x , 

3) A 3 =Y 3 = 2td t + xd x -2h8 h , 



4) A 4 = Y 4 = e~ u Ed E , 

5) A 5 = aY 1 +Y 2 =ad t + d x , 

6) A 6 = (3Y 3 + Y 4 = 2[3td t +f3xd x 



(4.34) 



2[3hd h + e- u Ed E , 



for nonzero constants a, (3. 



Proof. Let g 4 is the symmetry algebra of Eq. (1.1) with adjoint representation determined in Table 2 and 

Y = a 1 Y 1 +a 2 Y 2 + a 3 Y 3 + a 4 Y 4 (4.35) 
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is a nonzero vector field of g. We will simplify as many of the coefficients a» as possible through proper adjoint 
applications on Y. We follow our aim in the below easy cases. 



Case 1 At first, assume that 0,4 0. Scaling Y if necessary, we can consider 04 to be 1 and so follow the problem 
with 

Y = 01 Y 1 +a 2 Y 2 + a 3 Y 3 + Y 4 . (4.36) 
According to Table 2, if we act on Y by Ad(exp(^ a\ Y\)), the coefficient of Y\ can be vanished: 

Y' = a 2 Y 2 +a 3 Y 3 + Y 4 . (4.37) 
Then we apply Ad(exp(a 2 Y 2 )) on Y' to cancel the coefficient of Y 2 : 

Y" = a 3 Y 3 + Y 4 . (4.38) 

Case la If a 3 ^ then we can not simplify the coefficient of Y 3 to be either +1 or — 1. Thus any one-dimensional 
subalgebra generated by Y with a 3 , 04 ^ is equivalent to one generated by /3Y 3 + Y4 which introduce parts 6) 
of the theorem for constant (i ^ 0. 

Case lb For a 3 — we can see that each one-dimensional subalgebra generated by Y is equivalent to ( Y4 ) . 

Case 2 The remaining one-dimensional subalgebras are spanned by vector fields of the form Y with 04 = 0. 

Case 2a If a 3 ^ then by scaling Y, we can assume that a 3 = 1. Now by the action of Ad(exp(i ai Yi)) on Y, we 
can cancel the coefficient of Y\ : 

Y = a 2 Y 2 + Y 3 . (4.39) 

Then by applying Ad(exp(a2 Y 2 )) on Y the coefficient of Y 2 can be vanished and we tend to part 3) of the theorem. 
Case 2b Let a 3 — then Y is in the form 

Y = a 1 Y 1 +a 2 Y 2 . (4.40) 

Suppose that a 2 ^ then if necessary we can let it equal to 1. the simplest possible form of Y is equal to 

Y = ai Yi + Y 2 after taking a 2 = 1. 
Case 2b-l Let a\ be nonzero. In this case we can not make the coefficient of Y\ in Y more simpler and find 5) 

section of the theorem. ^ 
Case 2b-2 If a\ is zero then by scaling we can make the coefficient of Y equal to 1. Hence this case suggests part 

2). 

Case 2c Finally if in the latter case a 2 be zero, then no further simplification is possible and then Y is reduced to 
!)■ 

There is not any more possible case for studying and the proof is complete. □ 

The coefficients E, h of Eq. (1.1) resp. depend on the variables u, x. Therefore, we take their optimal system's 
projections on the space (a;, u, E, h). The nonzero in x— axis or u— axis projections of (4.34) are 



1) Z 1 = A 2 = A 5 =d x , 3) Z 3 = A* = e- u Ed El 

2) Z 2 = A 3 = xd x - 2hd h , 4) Z 4 = A 6 =/3xd x - 2/3hd h + e~ u Ed E - 

From paper 7 of [10] we conclude that 



(4.41) 



Proposition 5 Let g m := ( Yi : i = 1, • • • , m) be an m- dimensional algebra. Denote by A 1 (i = 1, • • • ,s, < s < 
m, s G N) an optimal system of one-dimensional subalgebras of g m and by Z % (i = 1, • • • , t, < t < s, t G N) the 
projections of A 1 , i.e., Z l =pr(A l ). If equations 

F = F(u), g = g(x), (4.42) 
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are invariant with respect to the optimal system Z % then the equation 

u t = (F(u) u x ) x + g(x) u, (4.43) 
admits the operators X 1 = projection of A 1 on (t,x,u). 

Proposition 6 Let Eq. (4-43) and the equation 

u t = (F(u)u x ) x +g(x)u, (4.44) 

be constructed according to Proposition 5 via optimal systems Z 1 and Z 1 resp. If the subalgebras spanned on the 
optimal systems Z % and Z l resp. are similar in Q m , then Eqs. (4-4-3) and (4-44) are equivalent with respect to the 
equivalence group G m generated by Q m . 

Now by applying Propositions 5 and 6 for the optimal system (4.41), we want to find all nonequivalent equations 
in the form of Eq. (1.1) admitting £ -extensions of the principal Lie algebra Qg , by one dimension, i.e, equations of 
the form (1.1) such that they admit, together with the one basic operators (1.2) of jji, also a second operator X^ 2 \ 
In each case which this extension occurs, we indicate the corresponding coefficients E, h and the additional operator 

We perform the algorithm passing from operators Z l (i = 1, • • • , 4) to E, h and X^ via the following example. 
Let consider the vector field 

K + "BE + '"~*SS' (445) 

then the characteristic equation corresponding to Z 4 is 
dx dh dE 



(4.46) 



(3x 2f3h e- u E' 
which determines invariants. Invariants can be taken in the following form 

h =u, I 2 = h 1 ' 2 x, h = x-^H 3 E. (4.47) 

In this case there are no invariant equations because the necessary condition for existence of invariant solutions (see 
[14], Section 19.3) is not satisfied, i.e., invariants (4.47) cannot be solved with respect to E and h since I3 is not an 
invariant function of I\ and I 2 to derive a function in term of u for E. 

Considering Z 2 we have the below characteristic equation 

dx dh 

-x = 27? (4 - 48) 
This equation suggest the following invariants 

h = u, h = h 1/2 x, I 3 = E. (4.49) 

From the invariancc equations we can write 

h = <t>{h), h = i>(h), (4.50) 

provided that <p(h) is a function in term of variable u and ip(Ii) an arbitrary function. The first condition occurs 
when (f){I\) = c for constant c. It results in the forms 

E = 4>(u), h ={^) 2 - ( 4 - 51 ) 
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Table 3 

The result of the classification 



N 


Z 


Invariant 


Equation 


Additional operator 


1 

2 


Z 1 

z 2 


u 
u 


u t = (4>(u) U x ) x + CU 
u t = {4>{u)u x ) x + (f) 2 w 


dx ' " at ~ dx 
It — + r — 

ZL at ^ x ax 



From Proposition 5 applied to the operator Z 4 we obtain the additional operator 



One can perform the algorithm for other Z l s of (4.41) similarly. The preliminary group classification of nonlinear 
fin equation (1.1) admitting an extension g 2 of the principal Lie algebra jji is listed in Table 3. 

5 Conclusion 

Projective analysis as a new symmetry property of equations ut = (E(u) u x ) x + h(x) u rather than previous results on 
this equation [17,18], is carried out exhaustively. Also, equivalence classification is given of the equation admitting an 
extension by one of the principal Lie algebra of the equation. The paper is one of few applications of a new algebraic 
approach to the problem of group classification: the method of preliminary group classification. Derived results are 
summarized in Table 3. 
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